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ABSTRACT 

We explore the hypothesis that the magnetic fields of neutron stars are of fossil 
origin. For parametrised models of the distribution of magnetic flux on the Main 
Sequence and of the birth spin period of the neutron stars, we calculate the expected 
properties of isolated radio pulsars in the Galaxy using as our starting point the initial 
mass function and star formation rate as a function of galacto-centric radius. We then 
use the 1374 MHz Parkes Multi-Beam Survey of isolated radio pulsars to constrain the 
parameters in our model and to deduce the required distribution of magnetic fields 
on the main sequence. We find agreement with observations for a model with a star 
formation rate that corresponds to a supernova rate of 2 per century in the Galaxy 
from stars with masses in the range 8 — 45M and predict 447,000 active pulsars in the 
Galaxy with luminosities greater than 0.19 mJy kpc 2 . The progenitor OB stars have 
a field distribution which peaks at ~ 46 G with ~ 8 percent of stars having fields in 
excess of 1,000 G. The higher field progenitors yield a population of 24 neutron stars 
with fields in excess of 10 14 G, periods ranging from 5 to 12 seconds, and ages of up 
to 100,000 years, which we identify as the dominant component of the magnetars. We 
also predict that high field neutron stars (log-B > 13.5) originate preferentially from 
higher mass progenitors and have a mean mass of 1.6Mq, which is significantly above 
the mean mass of 1.4M Q calculated for the overall population of radio pulsars. 
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1 INTRODUCTION 

A striking feature of the white dwarfs and the neutron stars 
is the wide range of magnetic field strengths seen in each 
of these groups. In isolated magnetic white dwarfs (about 
15 percent of the total white dwarf population), the field 
strengths are measured directly through the Zeeman effect 
and range from ~ 10 5 — 10 9 G. In isolated radio pulsars, 
the field estimates are in the range ~ 10 11 — 10 14 G and are 
based on the measured spin down rates and the assumption 
of dipole radiation. Observations in the X and 7-ray spec- 
tral regions have revealed the presence of a class of neutron 
stars with higher inferred fields, namely the Anomalous X- 
ray Pulsars (AXPs) and the Soft Gamma Repeaters (SGRs). 
These neutron stars are often referred to as "magnetars" and 
have estimated field strengths in the range ~ 10 14 - 10 15 G. 

Possible explanations for the magnetic fields in neutron 
stars are that the fields are either (i) essentially of fossil 
origin (e.g. Rudermann 1972) or (ii) generated by a connec- 
tive dynamo (Thompson & Duncan 1993). At the present 
time, this issue remains unresolved (see Mestel, 1999, for a 
comprehensive review). In the case of the white dwarfs, a 



direct link can be made between magnetism in the compact 
star phase and the Main Sequence via the well documented 
properties of the chemically peculiar Ap and Bp stars (Wick- 
ramasinghe & Ferrario 2005) and has been a major factor 
supporting (i) for the white dwarfs. 

Magnetism is not as well documented in the more mas- 
sive (^8M ) upper Main Sequence stars that evolve into 
neutron stars. The only known stars with a directly mea- 
sured magnetic field are 8 Orion C, with a dipolar magnetic 
field B p = 1,100 G (Donati et al. 2002), and HD191612 
with a field of B p = 1,500 G (Donati et al. 2005). Their 
mass is estimated around 40Mq , with the latter star slightly 
more evolved than the former. Interestingly, the magnetic 
fluxes of both these stars (1.1 x 10 27 G cm 2 for 9 Orion C 
and 7.5 x 10 27 G cm 2 for HD191612) are comparable to the 
flux of the highest field magnetar SGR 1806-20 (5.7 x 10 27 
G cm 2 ). There is also indirect evidence from studies of 
stellar winds from massive stars that magnetism may be 
widespread among OB stars (e.g. Wade 2001). These facts, 
taken together with the scaling that exists between the dipo- 
lar field strengths of the white dwarfs and the neutron stars 
(Figure 1), suggest that a prima facie case can be made for a 
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Figure 1. Black histogram (solid line): Field distribution of neu- 
tron stars. Red histogram (dashed line): Field distribution of High 
Field Magnetic White Dwarfs when their radii are shrunk to that 
of a typical neutron star radius of ~ 10 6 cm. 

fossil origin of fields also in the neutron stars. The fossil field 
hypothesis provides a natural explanation for the wide range 
of observed magnetic field strengths in a given class as be- 
ing the result of inhomogeneities in the interstellar magnetic 
field strength in star forming regions. 

Calculations of pre-Main Sequence evolution have 
shown that stars more massive than ~ 2M0 are likely to be- 
gin their Main Sequence phase with a primordial magnetic 
flux entrapped in radiative regions supporting the fossil hy- 
pothesis (Tout, Livio & Bonnell 1999; Moss 2003). There 
are no detailed studies on the interplay between fossil fields 
and subsequent stellar evolution. This problem has been dis- 
cussed by Tout, Wickramasinghe & Ferrario (2004, hereafter 
TWF). They argued that provided there are regions of the 
star that are radiative during subsequent evolution, a fossil 
Main Sequence poloidal magnetic flux could survive through 
to the compact star phase. 

In this paper, we investigate the consequences of the 
fossil field hypothesis for the origin of magnetic fields in neu- 
tron stars by carrying out population synthesis calculations 
for different assumptions on the distribution of the magnetic 
flux of massive ( ~ 8M©) Main Sequence stars and of the field 
dependence of the initial birth period of neutron stars. We 
model the observed properties of the population of isolated 
radio pulsars in the 1374 MHz Parkes Multi-Beam Survey 
(PMBS, Hobbs et al. 2004; Kramer et al. 2003; Morris et 
al. 2002; Manchester et al. 2001) to constrain our model 
parameters, and use these to deduce the required magnetic 
properties of the progenitor OB stars. 



2 THE MODEL 
2.1 The progenitors 

We assume that neutron stars are formed in the galactic 
disc from stars with 8 < M/Mq < 45 and that black holes 
form from larger mass stars. The situation is not as clear 
cut, from a theoretical point of view, because of the uncer- 
tainties in estimating the effects of fall-back (e.g. Heger et 
al. 2003). Nonetheless, there is now overwhelming observa- 
tional evidence that progenitors of neutron stars could be as 



massive as 40 — 50Mq (e.g. Gaensler et al. 2005; Muno et 
al. 2005). 

We assume an initial mass function if) (Mi) with a power 
law index Q3 = 2.7 (Kroupa 2002), and a constant star for- 
mation rate S(r g ) at a given galacto-centric radius r g . The 
radial dependence of S(r g ) is obtained from the data pre- 
sented by Boissier et al. (2003). The stars are distributed 
perpendicularly to the galactic plane (the z direction) expo- 
nentially using the scale heights as a function of initial mass 
given by Rana (1987) supplemented at the high mass end 
by the observed Wolf-Rayet distribution of Conti & Vacca 
(1990). These values range from a scale height of 100 pc for 
a 8.4M0 star to 45 pc for a 45M© star, which is in agreement 
with the more recent results of Reed (2000) . 

Unlike in the case of the white dwarfs, for the neu- 
tron stars the initial (Mi) to final mass (Mas) relationship 
is not well constrained by observations. Masses are known 
for a handful of neutron stars, particularly in binaries, and 
these are observed to cluster around Mns ~ 1.4Mq, close 
to the core collapse mass for an initial mass Mi ~ 8M0. In 
our calculations, we use the relationship obtained by Heger, 
Woosley & Spruit (2005). 

In the simple model of TWF, a star of mass Mi that 
collapses from the interstellar medium entraps a magnetic 
flux that is proportional to the ambient interstellar (primor- 
dial) magnetic field Bism and the square of its initial radius 
Ri. If we assume an initially nearly uniform density, it fol- 
lows that the magnetic flux <E> oc BismM^^ 3 . The variations 
expected in the interstellar magnetic field will then deter- 
mine the distribution of magnetic fluxes in the progenitor 
stars. We parametrise the magnetic flux distribution x(^) 
by assuming that all massive stars are magnetic, and that 
the magnetic flux $ is given as the sum of two Gaussians in 
the logarithm with dispersions CTiog*i and ffi g*2 with ap- 
propriate weightings. The mean values are assumed to have 
the form 

(log •)= log (*„.) + ! Jog (^) (1) 

where the subscript m takes the values 1 and 2 for the two 
Gaussians. The magnetic fluxes are in units of G cm 2 . 

2.2 Parametrisation of the birth properties of the 
neutron stars 

With the prescription of §2.1, we generate neutron stars with 
a range of masses. Our assumption that the magnetic flux is 
conserved during post main sequence evolution to the com- 
pact star phase allows us to calculate the birth surface mag- 
netic fields of these neutron stars from: 

Bns = * G (2) 
Tr(rtNs) 

where J?ns is the radius of the neutron star in cm given 
by the mass-radius relationship for neutron stars which de- 
pends on the equation of state. We adopt the values obtained 
with the unified SLy equation of state of Douchin & Haensel 
(2001). 

We follow the motions of these stars by integrating the 
equations of motion in the galactic potential of Kuijken & 
Gilmore (1989) assuming that they are born with a kick 
velocity that is independent of the progenitor mass and is 
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given by a Gaussian distribution (e.g. Hobbs et al. 2005) 
with velocity dispersion a v . 

Radio pulsars are generally observed not to have a brak- 
ing index of 3 (e.g. Cordes & Chernoff 1998). However, to 
minimise the number of free parameters we assume, in com- 
mon with many other theoretical studies, that the spin down 
of radio pulsars is by dipole radiation: 



described by a mean luminosity of the form 
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(3) 

where / is the moment of inertia as derived by Douchin & 
Haensel (2001) and P is the period in seconds. We have ex- 
plicitly assumed that the time scale for field decay is longer 
than the ages of the pulsars, although the issue of whether 
field decay is necessary to explain their detailed properties 
remains an open question (e.g. Gonthier et al. 2002). 

We may expect the birth period of the compact star to 
be determined by the magnetic flux of the progenitor star 
and its initial angular momentum (e.g. Brecher & Chan- 
mugam 1983; Narayan 1987). The complex interactions that 
occur between magnetic fields and rotation during stellar 
evolution and the mechanisms responsible for the trans- 
port of angular momentum are still poorly understood (e.g. 
Spruit & Phinney 1998; Watts & Andersson 2002; Fryer & 
Warren 2004; Heger et al. 2005; Ott et al. 2005). Magnetic 
fields must play a role in transporting angular momentum 
outwards from the stellar core to the stellar envelope be- 
cause otherwise all compact stars would be born at break 
up velocity, which is not the case (e.g. Migliazzo et al 2002 
and Kramer et al. 2003). Heger et al. (2005) have consid- 
ered this problem, but only in the weak field regime not 
directly relevant to the fossil field hypothesis. They report 
that the rotation periods they find are generally too small 
to explain pulsar observations. On the other hand, Ott et al. 
(2005) have proposed that a propeller mechanism operating 
on fall-back material may slow down the rotation of pulsars 
at birth. Interestingly, we know that there is strong direct 
evidence for higher field objects to be slower rotators in the 
magnetic white dwarfs and we speculate that this may be 
a general characteristic of stars with magnetic fields of fos- 
sil origin (see Ferrario & Wickramasinghe 2005 for scaling 
of rotation properties) and thus allow for this possibility in 
our modelling. 

We model the birth spin periods of neutron stars 
Po(Pns) by assuming they are distributed as a Gaussian 
n(Po) with a dispersion ap a about a mean Po(Bns) which 
depends linearly on the magnetic field 
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This is clearly a first approximation (e.g. in fall back mod- 
els, a field dependence oc B^J may be more appropriate). 
For different assumed values of a > 0, we constrain the pa- 
rameters P a , and crp by fitting the observations of isolated 
radio-pulsars in the PMBS taken from the ATNF catalogue 
(Manchester et al. 2005). 

To fit these observations, we also have to model the 
radio luminosity. Guided by the expected power law depen- 
dence on P and P, and following previous investigators, we 
have assumed that the luminosity L400 at 400 MHZ can be 



log(L 4 oo) = - log 



pS 



log L 



(•») 



with log Lo = 7.2. Here the luminosities are in units of mjy 
kpc 2 . We have modelled the spread around L400 using the 
dithering function of Narayan & Ostriker (1990). 

pi, (A) = 0.5A 2 exp(-A) (A > 0) 



where 



A = b log 



-J/400 
(£400 ) 



+ 
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and a and b are constants to be determined (Hartman et 
al. 1997). The luminosity at 400 MHZ is then scaled to the 
observed PMBS frequency of 1374 MHZ using a spectral 
index of -1.7 (Gonthier et al. 2002). 

In our studies, we have adopted the radio death line 
predicted for a multipolar field configuration near the stellar 
surface in a space-charge-limited flow model (Zhang, Hard- 
ing & Muslimov 2000), that is: 



log 



2 log 
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(8) 



2.3 The calculation of population characteristics 

We calculate the total number of neutron stars in the Galaxy 
with period up to P by evaluating the integral 

N{<P)= / / / / / ^(M l )x(M l ,$)n(P )r 9 5( 

JO JO J P 1 J PqJ 'S'minJ & 

dP 

x dMi d$ — dP dr„ dd> 
p 

Here, r g (in kpc) is the galacto-centric radial coordinate, and 
4> is the galactic longitude. P is calculated from equation [3] 
in appropriate units using the magnetic field Bj<ss(^{Mi)) 
derived from equation |2| The limits P01, P02, <&mm, $mai 
are chosen to adequately sample the Gaussian distributions. 

To compare our theoretical calculations with the PMBS 
isolated pulsars, we need to know the current location of 
each star that can contribute to the integral. To achieve 
this, we distribute the stars in velocity space according to 
the prescription of §2.2 with a v = 380 km s _1 and in the 2 
direction according to their progenitors masses as detailed 
in §2.1. Their spatial location at the current epoch is then 
determined by solving the equations of motion in the galac- 
tic potential. The dispersion and scattering measures DM 
and SM of the pulsars are calculated using the FORTRAN 
program of Cordes & Lazio (2002). Finally, the pulsar lumi- 
nosities are assigned according to equations to |7| 

Once all the intrinsic properties of our model pulsars 
are determined, these are fed through a filtering program 
kindly provided to us by Natasa Vranesevic (2005, private 
communication) that checks the pulsars for detectability by 
the PMBS. 

Finally, since pulsars radio emission is anisotropic, we 
use the Tauris & Manchester (1998) beaming model to ob- 
tain the correct fraction of detected pulsars: 



/(P) = 0.09 



log 



+ 0.03 



(10) 
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3 RESULTS AND DISCUSSION 

Our approach has been to focus on constraining the key 
parameters that we have introduced for describing the mag- 
netic flux distribution on the Main Sequence and the distri- 
bution of the initial birth period assuming that the remain- 
ing parameters of the model are within the ranges given by 
previous investigators. We have used as our basic observa- 
tional constraints the 1-D projections of the data comprising 
the number distributions in period P, magnetic field -Bjvs, 
period derivative P and radio luminosity L1374, noting that 
these distributions are not all independent. 

The best fit model to the observations of the PMBS 
pulsars was determined "by eye" after conducting hundreds 
of trials. We have found this method more reliable than the 
standard K-S statistic. Our results are shown in Figure 2. 

The initial period distribution is described by P a = 
0.22 s, a Pg = P /2 and a = 0.01. The magnetic flux 
distribution is described by two Gaussians of means and 
spreads log$i = 26.2, oi og 4>i = 1.0, and log $2 = 25.2, 
ciog $2 = 0.61 respectively, weighted in the ratio 1 : 5. The 
parameters for the luminosity model are a — 1.8 and 6 = 2.8. 

This model reproduces the observed total number of 
pulsars in the PMBS, once we exclude the millisecond pul- 
sars (P < 20 ms). The star formation rate is normalised to 
2.5M© pc -2 Gyr -1 at the galacto-centric radius of the Sun, 
which is within the uncertainty in the Boissier et al. (2003) 
relationship. This corresponds to a galactic supernova rate 
of 2 per century. 

In a broad sense, the distribution that we use for the 
magnetic flux has the strongest influence on the location 
of the peak in the magnetic field distribution of the pul- 
sars, while the mean period P a and the width of the pe- 
riod distribution o\p have a strong influence on the period 
distribution, although clearly period and field distributions 
are inter-linked through period evolution which depends on 
the magnetic field. We find that we cannot model both the 
field distribution and the period distribution at the same 
time if we assume that all pulsars are born at a period less 
than 50 ms. Such models give an excess of stars at low pe- 
riods which is not consistent with the results of Vranesevic 
et al. (2004) who find that up to 40 percent of pulsars are 
born with periods in the range 0.1-0.5 s. Regimbau & Freitas 
Pacheco (2001) also noted that most pulsars are born with 
P > 0.1 s and conduct their population synthesis studies us- 
ing a Gaussian with a mean birth period Po — 0.29 s and a 
spread of 0.1 s while Faucher-Giguere & Kaspi (2005) use a 
Gaussian with a mean birth period Po = 0.3 s and a spread 
of 0.15 s. All these results contradict the view that all pul- 
sars are born as fast rotators (Po < 0.1 s) which arises from 
observational selection effects, since young and fast neutron 
stars (e.g. the Crab pulsar) would be more easily detectable 
as radio pulsars than the more slowly rotating ones. We note 
that with our model parameters, the pulsar birth spin period 
increases with increasing magnetic field, but is almost inde- 
pendent of the field in the range log Bns(G) = 10—13 where 
the vast majority of isolated radio pulsars lie. The magnetic 
field dependence becomes important in the high-field radio 
pulsars and magnetar field range. 

The observations show that the field distribution is in- 
trinsically asymmetric with a low field tail that is more 
prominent than the high field tail. This general behaviour 
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Figure 3. Masses of neutron stars with rotational periods of up 
to 20 s in the two indicated field ranges. 

is reproduced by our assumed bi-Gaussian flux distribu- 
tion. Furthermore, the model predicts 26 pulsars with fields 
> 10 13 G, in good agreement with the PMBS which shows 
23 stars in the same field range. 

There have been many previous attempts at synthesis- 
ing the population of radio pulsars with different assump- 
tions on neutron star magnetic fields, field decay and initial 
periods. Some of these investigations have led to the con- 
clusion that the neutron star field distribution cannot be 
described by a single Gaussian (e.g. Gonthier et al. 2002), 
and we find a similar result for the magnetic fluxes of their 
progenitor stars. 

The initial-final mass radius relationship that we have 
assumed implies that higher mass progenitors will produce 
higher mass neutron stars. These stars will, on average, also 
have higher magnetic fluxes, and therefore tend to produce 
higher field neutron stars. However, the latter effect is coun- 
terbalanced by the strongly declining initial mass function 
which favours low mass star progenitors, which can also pro- 
duce high field neutron stars in the Gaussian tail of the 
magnetic flux distribution. In order to investigate the mass 
distribution of active radio pulsars in the Galaxy with pe- 
riods less than 20 s, we have divided them into two field 
groups (log Bjvs (<?) = 10-13.5) and (logBjvs(G) > 13.5). 
Our results are shown in Figure 3 for pulsars with 1374 
MHZ luminosities > 0.19 mjy kpc 2 . The mean neutron star 
masses in these groups are 1.45Mq and I.6M0 respectively, 
whilst the mean masses of their progenitors are 12M© and 
19M0 respectively. The stars in the highest field bin have a 
high mass tail with 54 percent of the stars extending beyond 
1.5Mq compared to 19 percent of stars in the same mass 
range in the lower field bin. We find that the total number of 
active pulsars with 1374 MHZ luminosities > 0.19 mjy kpc 2 
is 447,000. If we apply the beaming model of equation 1101 
we find that the number of active pulsars whose beams in- 
tercept Earth reduces to 48,000. 

The magnetars are neutron stars that are discovered as 
X-ray and/or 7-ray sources. They are currently confined to 
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Figure 2. Our model (blue) overlapped to the PMBS isolated radio pulsars data from the ATNF catalogue 
( http://www.atnf.csiro.au/research/pulsar/psrcat I. The errors are given by %/]V, where JV is the number of objects in each bin. 



the period range 5 — 12 s and are young objects (e.g. Gaensler 
et al. 2001). If one makes the standard assumption of dipole 
radiation with a braking index of 3, the magnetic fields are 
in the range log£>Ns(G) = 13.8 — 15.3. 

We have used our calculations to investigate if the high 
field tail of the neutron stars in our model can provide an 
explanation for the numbers of magnetars that are seen in 
the Galaxy. Here we use as our cutoff Bns = 10 14 G, since 
this value is just above the revised quantum critical field of 
Zhang & Harding (2000) for radio emission. If we isolate the 
total number of neutron stars with fields in excess of 10 14 G, 
and ages of less than 500,000 years, we find 146 such stars, 
none of which would be detectable as radio-pulsars. If we 
now limit the age to 100,000 years, our modelling yields 24 
magnetars in the Galaxy with periods from 5 to 12 s and 8 
magnetars with periods > 12 s, 5 of which have periods be- 
tween 12 and 15 seconds. Thus, the number declines rapidly 
with period outside the observed period range. We note in 
this context that some magnetars exhibit a large discrep- 
ancy between their real ages (as determined through their 
association with supernova remnants) and their characteris- 
tic ages (e.g. AXPs 1E2259+586 and 1E1845-0545, Gaensler 
et al. 2001). This discrepancy is readily explained through 
our assumption of a P — B relationship, which sees high 
field neutron stars being born preferentially at long periods 
(1 - 10 s). 

Below Bns = 10 14 G, there is an area of overlap where 
a couple of magnetars have fields and periods that are com- 
parable to those of the high field radio pulsars (HBRPs) 
(e.g. Kaspi & McLaughlin 2005). This suggests that there 
could exist lower field magnetars (< 10 14 G) that eventu- 
ally evolve into X-ray silent HBRPs of the type currently 
observed if their radio-beams intercept the Earth. We note 
that it should also be possible to discover "hybrid" young 
objects with _Bns < 10 14 G exhibiting both magnetar and 
radio pulsar characteristics. In our analysis, we have counted 
all these objects as radio-pulsars if they pass through our 
PMBS filtering program. 

Finally, we return to the recent evidence for high pro- 
genitor masses in the magnetars IE 1048-5937 (Gaensler et 
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Figure 4. Predicted magnetic field distribution of massive stars 
(8 — 45M©) on the Main Sequence. 



al. 2005) and CXO J164710.2-455216 (Muno et al. 2005). 
In the dynamo model, rapid spin (~ 3 ms) is required to 
generate super-strong fields and Heger et al. (2005) argue 
that it is only the cores of the more massive stars that have 
this property. This has been used as an argument in favour 
of the dynamo model, at least for the magnetars. We note 
that these findings also support our fossil hypothesis, since 
we predict that magnetars, and more generally high field 
neutron stars, originate preferentially from high mass pro- 
genitors. 

We conclude by noting that the fossil field hypothesis 
as formulated by us, and which does not allow for mag- 
netic flux loss in the post main-sequence evolution, requires 
a very specific distribution of magnetic fields for massive 
Main Sequence stars. This is shown in Figure 4 for an as- 
sumed dipolar field structure. This distribution has a peak 
at ~ 46 G and low and high field wings extending from 
~ 1 G to ~ 10, 000 G with 8 percent of stars having fields 
in excess of ~ 1, 000 G. The latter group are the progenitors 
of the highest field magnetars. The detailed shape of the 
Main Sequence distribution shown in Figure 4 is potentially 
an observable quantity, and a prediction of our fossil field 
model. 
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